Superconducting resonators as beam splitters for linear-optics quantum computation 
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A functioning quantum computer will be a ma- 
chine that builds up, in a programmable way, 
nonclassical correlations in a multipartite quan- 
tum system. Linear optics quantum computation 
(LOQC) is an approach for achieving this 

function that requires only simple, reliable lin- 
ear optical elements, namely beam splitters and 
phase shifters. Nonlinear optics is only required 
in the form of single-photon sources for state 
initialization, and detectors. However, the lat- 
ter remain difficult to achieve with high fidelity. 
A new setting for quantum optics has arisen in 
circuit quantum electrodynamics (cQED) using 
superconducting (SC) quantum devices [3], and 
opening up the way to LOQC using microwave, 
rather than visible photons. Much progress is 
being made in SC qubits [4] and cQED: high- 
fidelity Fock state generation [5] and qubit mea- 
surements provide single photon sources and de- 
tection. Here we show that the LOQC toolkit 
in cQED can be completed with high-fidelity (> 
99.92%) linear optical elements. 

We propose and analyze a technique for producing a 
beam-splitting quantum gate between two modes of a 
ring-resonator SC cavity. The cavity has two integrated 
superconducting quantum interference devices (SQUIDs) 
that are modulated by applying an external magnetic 
field. The gate is accomplished by applying a radio fre- 
quency pulse to one of the SQUIDs at the difference of 
the two mode frequencies. Departures from perfect beam 
splitting only arise from corrections to the rotating wave 
approximation; an exact calculation gives a fidelity of 
> 0.9992. 

We analyze in detail here a realization of a beam split- 
ter; the phase shifter represents a simpler case that can 
be studied in exactly the same way. Of course, a beam 
splitter is a well known device in microwave technology, 
where it is known as the directional coupler. But we need 
to achieve the functioning of a beam splitter in a different 
way, because the microwave photons in cQED are never 
"flying" ; we must find a way to accomplish the action 
of beam splitting in place between photonic modes that 
remain resident in SC resonators. 

In fact we show here that high-fidelity beam splitting 
can be produced in exactly the ring-resonator device an- 
alyzed by Ref. (see Fig. la). This structure has two 



nearly degenerate fundamental modes, even and odd with 
respect to the horizontal midline, with frequencies a;i.2 
respectively. The beam-splitting action will take place 
between photons in these two modes. This ring is in- 
terrupted by two SQUIDs as shown. In Q it is shown 
that the nonlinearity of these SQUIDs is enough to en- 
able nondemolition measurements; but this nonlinearity 
is far too small to make a practical traditional quantum 
gate such as a ePHASE. We can neglect this nonlinearity 
here, and we will exploit another control available in this 
device: the effective inductance of these SQUIDs can be 
controlled by the external magnetic fluxes shown. This 
inductance couples to the photonic modes by setting the 
reflection coefficient for standing waves around the ring. 
It is by modulation of the inductance of SQUIDi at the 
difference frequency uji —UJ2 that high-fidelity beam split- 
ting action is achieved 

We set up an accurate analysis of the photonic dynamics 
of our ring resonator by representing the rest of the ring 
connected to SQUIDi as a linear, passive structure with 
two-terminal impedance Z[(jj). Note that we will consider 
the flux bias of SQUID2 to be fixed, so that it functions 
as an inductor with, using the parameters of 6], Lq = 
0.19 nH. We neglect losses so that Z is purely imaginary. 
For device parameters from [6], Z{ljj) can be calculated 
analytically using standard two-port theory [s']. ln\Z{uj) 
in the frequency range of interest is shown in Fig. Ic. 
Since we include in Z the inductance Lq representing 
the average of the modulated SQUIDi, the two zeros of 
Z[(jj) correspond to the two mode frequencies uji^2 of the 
unmodulated ring resonator. 

This impedance is extremely well reproduced (dashed 
line) by that of the two-pole structure in our equivalent 
circuit shown in Fig. lb. This flt is obtained with param- 
eters Li = 5.8 nH, Ci = 86 fF, L2 = 7.7 nH, and C2 = 63 
fF. It is practical for the amplitude of the parametrically 
modulated inductance K{t) to be around 0.04 nH. We 
can apply our network graph theory Q to analyze the 
quantum behavior of these modes of the ring resonator. 
The classical time-dependent Hamiltonian that describes 
the equivalent circuit Fig. [T}) is 
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We choose the time dependent inductance to have the 
form Kit) = SK cos{uJdt). We will present calculations 
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FIG. 1: a) Schematic representation of the superconducting 
ring resonator, comprising two SG transmission line segments 
coupled by two SQUIDs that can be externally controlled 
by two different applied magnetic fluxes. AC modulation of 
SQUIDi will accomplish the action of beam splitting between 
the two fundamental modes of the ring resonator. The device 
can optionally be tuned by a stub on the midline as shown, 
b) Equivalent circuit, consisting of two parametrically coupled 
LG resonators. The modulation of the SQUIDi impedance is 
represented by the time-dependent inductor Kit'). The LC 
circuit parameters are chosen so that the impedance of the 
two parallel resonators matches the two-terminal impedance, 
looking down into the AB port, of the ring resonator, in a fre- 
quency band including the two fundamental modes, c) The 
AB-port impedance of the resonator and of the equivalent cir- 
cuit. The unmodulated fundamental mode frequencies lo\,2 
are given by the two zeros of this impedance. 



only for the case where 8K is switched on from zero at 
t = 0, then switched off again a.t t — t. We will see 
that even for such an unshaped, square modulation pulse, 
the desired quantum gate operation can be achieved with 
excellent fidelity. We find that this fidelity is insensitive 
to details of pulse shaping. 

We quantize Eq. ([1]) by imposing commutation rules be- 



tween canonically conjugate variables [^i,Qj 



ihSi. 



and express them, using creation and annihilation op- 
erators and a, as $i = ai\fh(ai + a|)/V2, and 
= -iVh{a, - al)/V2<7r, with ct,; = {U/a)^/^. As- 
suming K{t) <C ii,2, the Hamiltonian becomes (h — 1) 
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resonant harmonic frequencies are 
A = {LlC2/LlCif/\ and f{t) = 
f cos{ujdt). We will consider driv- 



ing at resonance at transition frequency ujd = Aw = 
ui — UJ2- For the ring resonator we will always be in 
the regime Acj ^ wi , • 

To study the beam splitting action created by the K{t) 
pulse, we will calculate the time evolution of an initial 
state comprising a single photon in one of the two modes 
(mode 1), i.e., IV-i) = al|0), with |0) = |0)i|0)2. We aim 
for the final "beam-split" state j-i/)/) = -^(oi — 
We consider always a doubly-rotating frame with rota- 
tion at frequency uji_2 for modes 1,2 resp. We will study 
a fidelity J-{t) = |(V'/|Z^(t)|V'i)P which indicates how 
close to ideal beam-splitting our operation is. Here, 
U{t) = Texp(— i J dt'H{t')) denotes the evolution op- 
erator generated by the Hamiltonian Eq. ([2]) ■ The use of 
more general gate fidelities would not affect our conclu- 
sions. 

We present three approaches to the calculation of T{t). 
First we perform a naive rotating wave approximation 
(RWA) in which only time independent terms in the ro- 
tating frame are retained in H{t). An elementary calcula- 
tion gives J-rwa(t) = (1 -l-sin(2/r))/2. Given the small- 
ness of Aw, a second calculation is much more accurate, 
which retains additional terms in the rotating frame that 
oscillate at frequencies Aw and 2 Aw. Then the part of 
the Hamiltonian that generates the beam-splitting gate 
is 
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The terms oscillating at frequency 2Aw produce Bloch- 
Siegert oscillations (BSO) 10 1. For / <C Aw these can 
be treated perturbatively [ll|. By considering virtual 
transitions via the first two sidebands shifted in energy 
by ±2 Aw, the fidelity at first order in //Aw is 
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Within both the RWA and the BSO approximations, for a 
pulse of duration t = 7r/4/ we attain = 1: beam split- 
ting is perfect at these levels of approximation. There is 
a real difference between these two points of view; in the 
naive RWA the Bloch vector in the |01)-|10) space un- 
dergoes simple circular motion, while the BSO manifest 
themselves as a nutational motion of the Bloch vector 
(Fig. 2). But in both cases the Bloch vector arrives ex- 
actly at the X-axis. Note that in both cases the Rabi 
oscillation frequency ilj^, = 2/ can be simply expressed 
as 
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For the parameters of our ring resonator, this gives a very 
convenient value of VIb/2ti k, 20MHz. 
Unfortunately, the actual time evolution of our Hamilto- 
nian does not give 100% fidelity for beam splitting. Para- 
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FIG. 2: To good approximation, tiie quantum evolution is 
confined to the Bloch sphere defined by the |01) and |10) 
states. Counterrotating terms in the Hamihonian at fre- 
quency 2Auj perturb the simple Bloch-sphere (red) evolution 
by terms of order //Aw. Bloch-Siegert oscillations (blue) 
cause nutation of the Bloch vector superimposed on regular 
precession in the rotating frame. The device parameters for 
the ring resonator are as given in the text. 



metric time-dependent modulation causes mode squeez- 
ing; in other words, photon number is not conserved, 
and our evolution does not remain confined to the |01)- 
|10) Bloch sphere. To quantify this effect, we must do a 
third calculation that goes beyond any rotating wave ap- 
proximation. With only a modest amount of numerical 
effort, it is feasible to do an essentially exact calculation 
of our gate operation. All that is required is a 4 x 4 
matrix calculation of the Heisenberg operators ai^2{t) 
and a\ 2 {t) . This calculation begins by using the 
canonically conjugate quadratures ^ = (51,(72,^1,^2)^ 
that are related to the original fluxes and charges by 
($1, $2, Qi, (52)"^ — VhDa^, with the diagonal matrix 
Da = diag(cri, (72, l/ci, l/o'2). Because the Hamilto- 
nian governing the evolution is quadratic, the Heisenberg 
equations of motion for the canonical quadratures ^ are 
linear, 
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where A and are real 2x2 matrices with il = 
diag(wi,a;2), An = A^, A22 = l/A^ and A12 = A21 = 1. 
The general solution can be expressed in terms of S(t) — 
Texp J dt'E{t'), where S{T) is a 4 X 4 real symplectic 
matrix that satisfies S'^{t)JS{t) = J, with the 4x4 real 
antisymmetric matrix J having a 2 x 2 block structure, 
with J12 = — J21 = 1 and Jn = J22 — 0. 
The action of the evolution operator U (t) on the canon- 
ical quadratures ^ results in a matrix multiplication 
^(t) = W{t)iU{t) = S{t)i that connects the Heisenberg 
and the Schrodinger representation. A similar relation 
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FIG. 3: Fidelity of the beam-splitting gate as a function of 
pulse duration, using a naive rotating wave approximation 
(RWA), a rotating wave approximation including some addi- 
tional slow-rotating terms perturbatively, capturing the oc- 
currence of Bloch-Siegert oscillations (BSO), and with an es- 
sentially exact numerical calculation. Due to its simple per- 
turbative nature, the BSO fidelity can (and does) exceed 1. 
The fidelity can be further optimized by carefully adjusting 
the resonator frequencies and thus the phase of the BSO. The 
device parameters for the ring resonator are as given in the 
text, except that stub tuning is used to raise the mode split- 
ting to Ao; 260 MHz. 



holds for the field operators and a\. By defining the 
vector a = (ai, a2, a|, 02)"^, the connection between the 
Heisenberg and the Schrodinger representation reads as 
a(t) = S'(c)(i)a, with S'(c)(t) = EtS'(t)E, where S is the 
simple unitary matrix of the basis change ^ = Sa. 

Here we show that the real symplectic matrix S{t) con- 
tains all the information needed to calculate the fidelity 
F{t) = |(i/'/|W(i)|V'i)P exactly. Any real symplectic ma- 
trix S e Sp(4,R) admits a singular value decomposi- 
tion in terms of real orthogonal symplectic matrices, S — 
SlDSI, with SlSa = 1, SlJSo, = J, for a = L,R, and 
D — diag(Ki,K2, 1/ki, I/K2), which is unique up to a re- 
ordering of the diagonal entries of D [l^l . This decompo- 
sition of S induces a corresponding decomposition of the 
evolution operator U{S): U[S) = U{Sl)U{D)U'^ [Sr). 
The elements Sa {a ~ L,R) have the general 2x2 block 

form [Sa]ll = [Sa]22 = Xa and [5'a]l2 = -[S'a]21 = Ya, 

with Xa, Ya real 2x2 matrices such that Ua = X^ — iYa 
is a unitary 2x2 matrix 
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It follows that the unitary evolution U{Sa) associ- 
ated with Sa conserves the number of photons and 
does not mix the creation and annihilation operators: 
{Sa)aiU{Sa) = J2k=i 2[^a]jkak [ill- On the other 
hand, the diagonal matrix D represents an active term 
that introduces squeezing in the two modes. In terms 
of mode operators U{D) can be expressed using inde- 



pendent squeezing operators Hi(p) = exp[p(af 
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1({D) ^ U{ki, K2) = ni(— ln«;i)n2(— lnK2)- Its action 
on the mode annihilation operators is simply given by 
Ill{p)ailli{p) = tti cosh/9 — a| sinh p. These facts allow us 
to write the fidelity as 



T = \{'ff\u{K,,^2m\' 



(7) 



with the state I^E* 



b,)\n)/V2- 



fi 



ib+ + 5+)|00)/V2 + (6r + 



b^\20) - 62^102) with bf ^ ± 1)/2kj 

and 6j = [Ur]ij{[Ul]ij + i[Ul]2j)* for j = 1,2. Eq. 
([7|) can now be evaluated in closed form, because the 
two-mode squeezed vacuum U{ki, K2)\0) has an analyti- 
cal expression in the Fock space via the known relation 
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with Hn{0) the Hcrmite polynomial at the origin, which 
is zero for odd n and -ff„(0) = (— l)"/^n!/(n/2)! for even 
n. Therefore, after a singular value deomposition of S 
the fidelity is directly obtained as a simple function of 
the quantities and 6i. 

Fig. 3 shows our three calculations of the beam-splitting 
fidelity as a function of pulse duration time r. The de- 
vice parameters for the ring resonator are as given above, 
except that by introducing an electrical stub for tun- 
ing as in Fig. la, the frequency of the even-symmetry 
mode uji/2tt is lowered to around 6.93 GHz, so that 
Auj/2tt sa 260 MHz. We see that the evolution approxi- 
mately follows the smooth Rabi oscillation predicted by 
the naive RWA, but that there are appreciable Bloch- 
Siegert oscillations superimposed on this. The perturba- 
tive BSO calculation in fact comes very close to the exact 
evolution for our parameters. The exact caclulation gives 
an extremely high value of the fidelity: Tmax > 0.9992. 
Squeezing is very small because the parametric modu- 
lation frequency is very slow compared with the mode 
frequencies (Aw <^ ^1,2)- For the ring resonator without 
the stub this ratio is even smaller, since then Aw « 64 
MHz; but this device is awkward to use since then the 
Bloch Siegert oscillations become very large and are no 
longer well described perturbatively [17 1. 

To summarize, we see that for a very straightforward 
ring-resonator geometry, almost ideal beam-splitting be- 
tween two cavity modes can be readily achieved. The 
calculated fidelity of 0.9992 is not realistic; other lim- 
its such as resonator loss and 1/f noise would come into 
play at this level for the present state of the art. But 
our result shows that there is no intrinsic limit to accom- 
plishing effective beam splitting by device modulation. 
Finally, we mention that an effective "delay line" is also 



readily implemented in this device; by transiently chang- 
ing the DC bias fluxes of the two SQUIDs, in either an 
even or an odd fashion, either mode may be subjected to 
any desired phase shift. There will be even less intrinsic 
limitation on the fidelity of these operations. Thus, we 
see that the toolkit for linear optics quantum comput- 
ing is readily completable in superconducting microwave 
circuits. 
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